A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations.

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, genome sequencing (GS), RNA-seq, and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined GS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints matches the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2 Mb region on chromosome 9 and a SINE element insertion at the more distal breakpoint. Interestingly, this genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by both RNA-seq and Sanger sequencing on blood samples from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p24 segregating with a familial congenital heart defect, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

Structural variants in the Epb41l4a locus: TAD disruption and Nrep gene misregulation as hypothetical drivers of neurodevelopmental outcomes.

Structural variations are a pervasive feature of human genomes, and there is growing recognition of their role in disease development through their impact on spatial chromatin architecture. This understanding has led us to investigate the clinical significance of CNVs in noncoding regions that influence TAD structures. In this study, we focused on the Epb41l4a locus, which contains a highly conserved TAD boundary present in both human chromosome 5 and mouse chromosome 18, and its association with neurodevelopmental phenotypes. Analysis of human data from the DECIPHER database indicates that CNVs within this locus, including both deletions and duplications, are often observed alongside neurological abnormalities, such as dyslexia and intellectual disability, although there is not enough evidence of a direct correlation or causative relationship. To investigate these possible associations, we generated mouse models with deletion and inversion mutations at this locus and carried out RNA-seq analysis to elucidate gene expression changes. We found that modifications in the Epb41l4a TAD boundary led to dysregulation of the Nrep gene, which plays a crucial role in nervous system development. These findings underscore the potential pathogenicity of these CNVs and highlight the crucial role of spatial genome architecture in gene expression regulation.

Cell type-specific interpretation of noncoding variants using deep learning-based methods.

Interpretation of noncoding genomic variants is one of the most important challenges in human genetics. Machine learning methods have emerged recently as a powerful tool to solve this problem. State-of-the-art approaches allow prediction of transcriptional and epigenetic effects caused by noncoding mutations. However, these approaches require specific experimental data for training and cannot generalize across cell types where required features were not experimentally measured. We show here that available epigenetic characteristics of human cell types are extremely sparse, limiting those approaches that rely on specific epigenetic input. We propose a new neural network architecture, DeepCT, which can learn complex interconnections of epigenetic features and infer unmeasured data from any available input. Furthermore, we show that DeepCT can learn cell type-specific properties, build biologically meaningful vector representations of cell types, and utilize these representations to generate cell type-specific predictions of the effects of noncoding variations in the human genome.

Lampbrush chromosome studies in the post-genomic era.

Extraordinary extended lampbrush chromosomes with thousands of transcription loops are favorable objects in chromosome biology. Chromosomes become lampbrushy due to unusually high rate of transcription during oogenesis. However, until recently, the information on the spectrum of transcribed sequences as well as genomic context of individual chromomeres was mainly limited to tandemly repetitive elements. Here we briefly outline novel findings and future directions in lampbrush chromosome studies in the post-genomic era. We emphasize the fruitfulness of combining genome-wide approaches with microscopy imaging techniques using lampbrush chromosomes as a remarkable model object. We believe that new data on the spectrum of sequences transcribed on the lateral loops of lampbrush chromosomes and their structural organization push the boundaries in the discussion of their biological role. Also see the video abstract here: https://youtu.be/zexoHfzX9rM.

Function and Evolution of the Loop Extrusion Machinery in Animals.

Structural maintenance of chromosomes (SMC) complexes are essential proteins found in genomes of all cellular organisms. Essential functions of these proteins, such as mitotic chromosome formation and sister chromatid cohesion, were discovered a long time ago. Recent advances in chromatin biology showed that SMC proteins are involved in many other genomic processes, acting as active motors extruding DNA, which leads to the formation of chromatin loops. Some loops formed by SMC proteins are highly cell type and developmental stage specific, such as SMC-mediated DNA loops required for VDJ recombination in B-cell progenitors, or dosage compensation in Caenorhabditis elegans and X-chromosome inactivation in mice. In this review, we focus on the extrusion-based mechanisms that are common for multiple cell types and species. We will first describe an anatomy of SMC complexes and their accessory proteins. Next, we provide biochemical details of the extrusion process. We follow this by the sections describing the role of SMC complexes in gene regulation, DNA repair, and chromatin topology.

A Cre-LoxP-based approach for combinatorial chromosome rearrangements in human HAP1 cells.

Alterations of human karyotype caused by chromosomal rearrangements are often associated with considerable phenotypic effects. Studying molecular mechanisms underlying these effects requires an efficient and scalable experimental model. Here, we propose a Cre-LoxP-based approach for the generation of combinatorial diversity of chromosomal rearrangements. We demonstrate that using the developed system, both intra- and inter-chromosomal rearrangements can be induced in the human haploid HAP1 cells, although the latter is significantly less effective. The obtained genetically modified HAP1 cell line can be used to dissect genomic effects associated with intra-chromosomal structural variations.

Annotation of uORFs in the OMIM genes allows to reveal pathogenic variants in 5'UTRs.

An increasing number of studies emphasize the role of non-coding variants in the development of hereditary diseases. However, the interpretation of such variants in clinical genetic testing still remains a critical challenge due to poor knowledge of their pathogenicity mechanisms. It was previously shown that variants in 5'-untranslated regions (5'UTRs) can lead to hereditary diseases due to disruption of upstream open reading frames (uORFs). Here, we performed a manual annotation of upstream translation initiation sites (TISs) in human disease-associated genes from the OMIM database and revealed ∼4.7 thousand of TISs related to uORFs. We compared our TISs with the previous studies and provided a list of 'high confidence' uORFs. Using a luciferase assay, we experimentally validated the translation of uORFs in the ETFDH, PAX9, MAST1, HTT, TTN,GLI2 and COL2A1 genes, as well as existence of N-terminal CDS extension in the ZIC2 gene. Besides, we created a tool to annotate the effects of genetic variants located in uORFs. We revealed the variants from the HGMD and ClinVar databases that disrupt uORFs and thereby could lead to Mendelian disorders. We also showed that the distribution of uORFs-affecting variants differs between pathogenic and population variants. Finally, drawing on manually curated data, we developed a machine-learning algorithm that allows us to predict the TISs in other human genes.

A germline chimeric KANK1-DMRT1 transcript derived from a complex structural variant is associated with a congenital heart defect segregating across five generations.

Structural variants (SVs) pose a challenge to detect and interpret, but their study provides novel biological insights and molecular diagnosis underlying rare diseases. The aim of this study was to resolve a 9p24 rearrangement segregating in a family through five generations with a congenital heart defect (congenital pulmonary and aortic valvular stenosis, and pulmonary artery stenosis), by applying a combined genomic analysis. The analysis involved multiple techniques, including karyotype, chromosomal microarray analysis (CMA), FISH, whole-genome sequencing (WGS), RNA-seq and optical genome mapping (OGM). A complex 9p24 SV was hinted at by CMA results, showing three interspersed duplicated segments. Combined WGS and OGM analyses revealed that the 9p24 duplications constitute a complex SV, on which a set of breakpoints match the boundaries of the CMA duplicated sequences. The proposed structure for this complex rearrangement implies three duplications associated with an inversion of ~ 2Mb region on chromosome 9 with a SINE element insertion at the more distal breakpoint. Interestingly, this hypothesized genomic structure of rearrangement forms a chimeric transcript of the KANK1/DMRT1 loci, which was confirmed by RNA-seq on blood from 9p24 rearrangement carriers. Altogether with breakpoint amplification and FISH analysis, this combined approach allowed a deep characterization of this complex rearrangement. Although the genotype-phenotype correlation remains elusive from the molecular mechanism point of view, this study identified a large genomic rearrangement at 9p segregating with a familial congenital clinical trait, revealing a genetic biomarker that was successfully applied for embryo selection, changing the reproductive perspective of affected individuals.

Multilevel view on chromatin architecture alterations in cancer.

Chromosomes inside the nucleus are not located in the form of linear molecules. Instead, there is a complex multilevel genome folding that includes nucleosomes packaging, formation of chromatin loops, domains, compartments, and finally, chromosomal territories. Proper spatial organization play an essential role for the correct functioning of the genome, and is therefore dynamically changed during development or disease. Here we discuss how the organization of the cancer cell genome differs from the healthy genome at various levels. A better understanding of how malignization affects genome organization and long-range gene regulation will help to reveal the molecular mechanisms underlying cancer development and evolution.

Method for the Isolation of "RNA-seq-Quality" RNA from Human Intervertebral Discs after Mortar and Pestle Homogenization.

The problem of isolating high-quality total RNA from intervertebral discs has no recognized solution yet. This is due to the extremely low content of live cells in the samples and the voluminous intercellular matrix. A variety of published protocols focused on isolating RNA from articular cartilage have recommended the use of expensive equipment, enzymatic matrix cleavage, or cell culture. In our study, we used a combination of the traditional QIAzol protocol (Qiagen, Germany) and RNEasy column purification (Qiagen, Germany) to obtain high-quality RNA from post-surgical intervertebral disc fragments. Only a mortar and a pestle were used for grinding, making our method particularly accessible. The isolated RNA with a RIN of ~7 is suitable for studying the expression profile of chondrocytes in situ. RNA-seq analysis of three samples demonstrated cell type ratios to be mostly relevant to intervertebral disc tissues, with over 70% of the chondrocytes of the three subtypes having an admixture of blood-related cells.

Single-Cell DNA Methylation Analysis of Chicken Lampbrush Chromosomes.

DNA methylation is an essential epigenetic regulation mechanism implicated in transcription and replication control, developmental reprogramming, retroelements silencing and other genomic processes. During mammalian development, a specific DNA methylation pattern should be established in germ cells to allow embryonic development. Less is known about germ cell DNA methylation in other species. To close this gap, we performed a single-cell methylome analysis of chicken diplotene oocytes. We comprehensively characterized methylation patterns in these cells, obtained methylation-based chicken genome segmentation and identified oocyte-specific methylated gene promoters. Our data show that despite the formation of specific transcriptionally hyperactive genome architecture in chicken diplotene oocytes, methylation patterns in these cells closely resemble genomic distribution observed in somatic tissues.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects.

Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

Derivation of Ringed Seal (Phoca hispida) Induced Multipotent Stem Cells.

Induced pluripotent stem (iPS) cells have been produced just for a few species among order Carnivora: snow leopard, Bengal tiger, serval, jaguar, cat, dog, ferret, and American mink. We applied the iPS cell derivation protocol to the ringed seal (Phoca hispida) fibroblasts. The resulting cell line had the expression of pluripotency marker gene Rex1. Differentiation in embryoid body-like structures allowed us to register expression of AFP, endoderm marker, and Cdx2, trophectoderm marker, but not neuronal (ectoderm) markers. The cells readily differentiated into adipocytes and osteocytes, mesoderm cell types of origin. Transcriptome analysis allowed us to conclude that the cell line does not resemble human pluripotent cells, and, therefore, most probably is not pluripotent. Thus, we produced ringed seal multipotent stem cell line capable of differentiation into adipocytes and osteocytes.

Differential DNA Methylation of the IMMP2L Gene in Families with Maternally Inherited 7q31.1 Microdeletions is Associated with Intellectual Disability and Developmental Delay.

Most copy number variations (CNVs) in the human genome display incomplete penetrance with unknown underlying mechanisms. One such mechanism may be epigenetic modification, particularly DNA methylation. The IMMP2L gene is located in a critical region for autism susceptibility on chromosome 7q (AUTS1). The level of DNA methylation was assessed by bisulfite sequencing of 87 CpG sites in the IMMP2L gene in 3 families with maternally inherited 7q31.1 microdeletions affecting the IMMP2L gene alone. Bisulfite sequencing revealed comparable levels of DNA methylation in the probands, healthy siblings without microdeletions, and their fathers. In contrast, a reduced DNA methylation index and increased IMMP2L expression were observed in lymphocytes from the healthy mothers compared with the probands. A number of genes were upregulated in the healthy mothers compared to controls and downregulated in probands compared to mothers. These genes were enriched in components of the ribosome and electron transport chain, as well as oxidative phosphorylation and various degenerative conditions. Differential expression in probands and mothers with IMMP2L deletions relative to controls may be due to compensatory processes in healthy mothers with IMMP2L deletions and disturbances of these processes in probands with intellectual disability. The results suggest a possible partial compensation for IMMP2L gene haploinsufficiency in healthy mothers with the 7q31.1 microdeletion by reducing the DNA methylation level. Differential DNA methylation of intragenic CpG sites may affect the phenotypic manifestation of CNVs and explain the incomplete penetrance of chromosomal microdeletions.

A cookbook for DNase Hi-C.

BACKGROUND: The Hi-C technique is widely employed to study the 3-dimensional chromatin architecture and to assemble genomes. The conventional in situ Hi-C protocol employs restriction enzymes to digest chromatin, which results in nonuniform genomic coverage. Using sequence-agnostic restriction enzymes, such as DNAse I, could help to overcome this limitation. RESULTS: In this study, we compare different DNAse Hi-C protocols and identify the critical steps that significantly affect the efficiency of the protocol. In particular, we show that the SDS quenching strategy strongly affects subsequent chromatin digestion. The presence of biotinylated oligonucleotide adapters may lead to ligase reaction by-products, which can be avoided by rational design of the adapter sequences. Moreover, the use of nucleotide-exchange enzymes for biotin fill-in enables simultaneous labelling and repair of DNA ends, similar to the conventional Hi-C protocol. These improvements simplify the protocol, making it less expensive and time-consuming. CONCLUSIONS: We propose a new robust protocol for the preparation of DNAse Hi-C libraries from cultured human cells and blood samples supplemented with experimental controls and computational tools for the evaluation of library quality.

Erythrocytes 3D genome organization in vertebrates.

Generation of mature red blood cells, consisting mainly of hemoglobin, is a remarkable example of coordinated action of various signaling networks. Chromatin condensation is an essential step for terminal erythroid differentiation and subsequent nuclear expulsion in mammals. Here, we profiled 3D genome organization in the blood cells from ten species belonging to different vertebrate classes. Our analysis of contact maps revealed a striking absence of such 3D interaction patterns as loops or TADs in blood cells of all analyzed representatives. We also detect large-scale chromatin rearrangements in blood cells from mammals, birds, reptiles and amphibians: their contact maps display strong second diagonal pattern, representing an increased frequency of long-range contacts, unrelated to TADs or compartments. This pattern is completely atypical for interphase chromosome structure. We confirm that these principles of genome organization are conservative in vertebrate erythroid cells.

The Mutation Spectrum of Maturity Onset Diabetes of the Young (MODY)-Associated Genes among Western Siberia Patients.

Maturity onset diabetes of the young (MODY) is a congenital form of diabetes characterized by onset at a young age and a primary defect in pancreatic-ß-cell function. Currently, 14 subtypes of MODY are known, and each is associated with mutations in a specific gene: HNF4A, GCK, HNF1A, PDX1, HNF1B, NEUROD1, KLF11, CEL, PAX4, INS, BLK, KCNJ11, ABCC8, and APPL1. The most common subtypes of MODY are associated with mutations in the genes GCK, HNF1A, HNF4A, and HNF1B. Among them, up to 70% of cases are caused by mutations in GCK and HNF1A. Here, an analysis of 14 MODY genes was performed in 178 patients with a MODY phenotype in Western Siberia. Multiplex ligation-dependent probe amplification analysis of DNA samples from 50 randomly selected patients without detectable mutations did not reveal large rearrangements in the MODY genes. In 38 patients (37% males) among the 178 subjects, mutations were identified in HNF4A, GCK, HNF1A, and ABCC8. We identified novel potentially causative mutations p.Lys142*, Leu146Val, Ala173Glnfs*30, Val181Asp, Gly261Ala, IVS7 c.864 -1G>T, Cys371*, and Glu443Lys in GCK and Ser6Arg, IVS 2 c.526 +1 G>T, IVS3 c.713 +2 T>A, and Arg238Lys in HNF1A.

Postsynthetic On-Column 2' Functionalization of RNA by Convenient Versatile Method.

We report a universal straightforward strategy for the chemical synthesis of modified oligoribonucleotides containing functional groups of different structures at the 2' position of ribose. The on-column synthetic concept is based on the incorporation of two types of commercial nucleotide phosphoramidites containing orthogonal 2'-O-protecting groups, namely 2'-O-thiomorpholine-carbothioate (TC, as "permanent") and 2'-O-tert-butyl(dimethyl)silyl (tBDMS, as "temporary"), to RNA during solid-phase synthesis. Subsequently, the support-bound RNA undergoes selective deprotection and follows postsynthetic 2' functionalization of the naked hydroxyl group. This convenient method to tailor RNA, utilizing the advantages of solid phase approaches, gives an opportunity to introduce site-specifically a wide range of linkers and functional groups. By this strategy, a series of RNAs containing diverse 2' functionalities were synthesized and studied with respect to their physicochemical properties.

DNA barcoding reveals that injected transgenes are predominantly processed by homologous recombination in mouse zygote.

Mechanisms that ensure repair of double-strand DNA breaks (DSBs) are instrumental in the integration of foreign DNA into the genome of transgenic organisms. After pronuclear microinjection, exogenous DNA is usually found as a concatemer comprising multiple co-integrated transgene copies. Here, we investigated the contribution of various DSB repair pathways to the concatemer formation. We injected mouse zygotes with a pool of linear DNA molecules carrying unique barcodes at both ends and obtained 10 transgenic embryos with 1-300 transgene copies. Sequencing the barcodes allowed us to assign relative positions to the copies in concatemers and detect recombination events that occurred during integration. Cumulative analysis of approximately 1,000 integrated copies reveals that over 80% of them underwent recombination when their linear ends were processed by synthesis-dependent strand annealing (SDSA) or double-strand break repair (DSBR). We also observed evidence of double Holliday junction (dHJ) formation and crossing over during the concatemer formations. Sequencing indels at the junctions between copies shows that at least 10% of DNA molecules introduced into the zygotes are ligated by non-homologous end joining (NHEJ). Our barcoding approach, verified with Pacific Biosciences Single Molecule Real-Time (SMRT) long-range sequencing, documents high activity of homologous recombination after DNA microinjection.